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To explore the roles of dynamical model and statistical model in the description of multifrag¬ 
mentation in heavy-ion collisions at intermediate energies, the fragments charge distribntions of 
^®^An-|-^®^An at 35 MeV/u are analyzed by using the hybrid model of improved quantum molecular 
dynamics (ImQMD) model plus the statistical model GEMINI. We find that, the ImQMD model can 
well describe the charge distributions of fragments produced in central and semi-central collisions. 

But for the peripheral collisions of And-An at 35 MeV/u, the statistical model is required at the 
end of the ImQMD simulations for the better description of the charge distribution of fragments. 

By using the hybrid model of ImQMD+GEMINI, the fragment charge distribution of And-An at 35 
MeV/u can be reproduced reasonably well. The time evolution of the excitation energies of primary 
fragments is simultaneously investigated. 

PACS numbers: 24.10.Cn, 24.10.Pa, 25.70.Mn 


I. INTRODUCTION 

The heavy-ion nuclear reactions at intermediate and 
high energies provide a unique approach to study the nu¬ 
clear equation of state at high densities in terrestrial lab¬ 
oratories, and therefore attracted much attentions both 
in the experimental and theoretical fields [1-8]. Because 
of the very short time scale in heavy-ion collisions, it 
is difficult to record the dynamical process of reactions 
by the experimental equipment. What the experiments 
can provide are the phase space information of various 
particles produced in the reactions. To know what hap¬ 
pened during the reactions, especially before the sys¬ 
tem equilibrium, reliable dynamical approaches are cru¬ 
cial to extract physics from heavy-ion collision experi¬ 
ments. Many transport theories have been developed to 
explain the relevant experimental data. The models can 
be classified into two main categories: one-body theory 
and n-body theory. The popular one-body theories in¬ 
clude the time dependent Hartree-Fock [9-11] and the 
Boltzmann-Vlasov approach[12-14], e.g. IBUU, RBUU, 
pBUU, RVUU, GiBUU, et. al, in BV family. One of the 
popular n-body theories is the molecular dynamics (MD) 
approach, e.g. QMD, AMD, ImQMD, IQMD, UrQMD, 
CoMD, et. al, in MD family. With many efforts on the 
study of heavy-ion reactions, more and more observables 
are obtained and need to be explained by the transport 
model, the reliability of transport models faces a great 
challenge. One of the most challenging tasks for the n- 
body theory in nuclear physics is the description of the 
multifragmentation of heavy nuclei. The investigation of 
multifragmentation [15-20] is important for understand¬ 
ing the reaction mechanism in heavy-ion collisions. One- 
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body theories are not proper tools for the investigation of 
many-body correlations, because the formation of frag¬ 
ments is beyond the scope of these models. Another kind 
of theory is the statistical approaches, i.e. statistical or 
thermodynamical models. Such as canonical thermody¬ 
namical model (CTM) [18], statistical multifragmenta¬ 
tion model (SMM) [19], GEMINI model [21, 22], GEM2 
[23, 24] model and HIVAP code [25]. The statistical ap¬ 
proach is more appropriate to study the nuclei in equilib¬ 
rium state. So the statistical approach is used to describe 
the spectator in peripheral heavy-ion collisions. 

The n-body theory can well describe the multifrag¬ 
mentation process in central and semi-central heavy-ion 
reactions at intermediate and high energies, in which the 
large density fluctuations caused by compression and ex¬ 
pansion and high excitation energies of the system play 
a dominant role for the formation of fragments. How¬ 
ever, the n-body theory fails to self-consistently describe 
the peripheral heavy-ion collisions and spallation reac¬ 
tions with heavy nuclei because in this kind of reactions 
the density fluctuations and excitation energies are rela¬ 
tively small and the fission of composite system or heavy 
fragments could not be neglected. So the statistical anal¬ 
ysis at the end of dynamical simulations is necessary and 
important, since the transport models only contain classi¬ 
cal correlations, which is insufficient to correctly describe 
evaporation (where realistic density of states are needed) 
or fission (where quantum fluctuations are essential). On 
the other hand, low energy reactions require a simulation 
with very long time. In practice it is a massive mission 
to simulate the reaction of very heavy system with very 
large time scale. The hybrid model by dynamical plus 
statistical model is a good compromise to study this kind 
of subject. 

In this paper, we use the improved quantum molecu¬ 
lar dynamics (ImQMD) model plus the statistical model 
GEMINI to describe the multifragmentation process in 
heavy ion collisions. The paper is organized as follows. 
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In sec. II, we briefly introduce the model we adopted. 
In sec. Ill, we present some calculations about the ex¬ 
citation energy and charge distribution of fragments in 
the reactions of ^®^Au-|-^®^Au at incident energies of 
Ec.m. =35 MeV/u. Finally a brief summary is given in 
Sec. IV. 


II. THEORETICAL APPROACHES 


In the ImQMD model [26], as in the original QMD 
model [27], each nucleon is represented by a coherent 
state of a Gaussian wave packet. 
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where and Pi are the centers of the ith wave packet in 
the coordinate and momentum space, respectively. The 
one-body phase space distribution function is obtained 
by the Wigner transform of the wave function. The time 
evolution of and Pi for each nucleon is governed by 
Hamiltonian equations of motion 
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here, the kinetic energy T = ^ C/coui is the 
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Coulomb energy, and the local potential energy C/ioc = 
/ Vioc[p{r)]dr. Vioc is the nuclear potential energy den¬ 
sity functional that is obtained by the effective Skyrme 
interaction, which reads 
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In this work, we adopt the parameter set IQS [28] (see Ta¬ 
ble I) , which has been proposed and tested for describing 
the heavy-ion fusion reactions and the multifragmenta¬ 
tion process, in the previous works [28, 29]. 

At the end of the ImQMD calculations, clusters are 
recognized by a minimum spanning tree (MST) algo¬ 
rithm [27] widely used in the QMD calculations. In this 
method, the nucleons with relative momenta smaller than 
Pc and relative distances smaller than Rc are coalesced 
into the same cluster. In this work, Rc = 3.5 fm and 
Pc = 300 MeV/c are adopted. Then the total energy of 
each excited cluster is calculated in its rest frame and its 
excitation energy are obtained by subtracting the corre¬ 
sponding ground state energy from the total energy of the 
excited cluster. The information of excited cluster is in¬ 
put into the statistical decay model GEMINI to perform 
statistical decay stage calculations. 


The GEMINI is based on the well-known sequential- 
binary-decay picture that the individual compound nu¬ 
clei decay through sequential binary decays of all possi¬ 
ble modes, from emission of nucleons and light particles 
through asymmetric to symmetric fission as well as the 
7 emission, until the resulting products are unable to 
undergo any further decay. Fortran95 version (released 
on 2-May-2005) of GEMINI is used in this work. A 
very detailed description of the GEMINI can be found 
in Ref. [21, 22] and references therein. The GEMINI pa¬ 
rameters are chosen as following: The level density is 
taken as Grimes case B modified form {aden_type=—23)^ 
all asymmetric divisions is considered in fission mode 
{imf-option=2), the particle with Z < 5 is treated in 
light particle evaporation (Z_imf_min=5). The other pa¬ 
rameters are taken as the default value given by example 
in GEMINI document. 


HI. RESULTS AND DISCUSSION 

As the most important input of GEMINI model, the 
excitation energy determines the decay process of the pri¬ 
mary cluster formed in the dynamical stage. So the cal¬ 
culation of excitation energy should be checked carefully. 
We firstly prepare the two initial nuclei with their prop¬ 
erties of the ground states being well described. The time 
evolutions of the binding energies and nuclear radii for a 
number of nuclei have been checked simultaneously. Only 
those individual nucleus which can remain stable for sev¬ 
eral thousands fm/c are taken to be good initial nuclei, 
and then are applied in the simulation of the reaction 
process. 

Figure 1 presents the excitation energy distribution for 
the head-on collisions of ^°Ga-|-^°Ca {Ec.m. = 80 MeV), 
i6o-^46Ti (^Ec.m. = 38 MeV), and {Ec.m. = 50 

MeV), at 200 fm/c and 500 fm/c, respectively. To check 
the procedure of the calculation of excitation energy in 
the ImQMD model, the corresponding excitation energy 
Eq = Ec.m.+Qgg of the compound nuclei in heavy-ion fu¬ 
sion reaction is also presented for comparison. Here Qgg 
denotes the Q value of the fusion reaction from ground 
state to ground state. Qgg = —14.3, 12.3, and —3.9 MeV 
for ^°Ca-|-^°Ca, ^®0-|-''^®Ti and ^®0-|-®^Zr, respectively. 
The corresponding values of Eq are 65.7, 50.3 and 46.1 
MeV for the three reaction systems (see the dashed lines), 
respectively. One can see that at t =200 fm/c, the pro¬ 
jectile and target are well separated, so the excitation 
energies of target nuclei are distributed around zero. At 
t =500 fm/c, the compound nuclei are generally formed 
and the peaks of the excitation energies distribution lo¬ 
cate at around the corresponding Eq. These tests indi¬ 
cate that the calculation results of excitation energies of 
primary fragments are reasonable. 

The reactions of ^®^Au-|-^®^Au at incident energies of 
Ec.m. = 35 MeV/u are simulated with the ImQMD model 
to investigate the multifragmentation behavior of heavy- 
ion collision. Firstly, we check the time evolutions of exci- 
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TABLE I: Parameter set IQS used in the ImQMD calculations. 


a d 
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(MeV fm^) 
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-207 138 
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18.0 

14.0 

5/3 

32.0 0.08 

0.165 

0.94 

0.018 
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FIG. 1: (Color online) Excitation energy distribution of the largest fragments for head-on collisions of ^*’Ca-|-^°Ca {Ec.m. = 80 
MeV) (left panel), ^®0+"^®Ti (Ec.m. = 38 MeV) (middle panel), and ^®0-|-®^Zr {Ec.m. = 50 MeV) (right panel), at 200 fm/c 
(top panel) and 500 fm/c (bottom panel), respectively. The contour plot of nuclear density distribution is also presented. 


tation energies per nucleon of primary fragments formed 
in reactions with impact parameters b=l, 6, and 10 fm, 
respectively, as presented in Fig. 2. One can get a di¬ 
rect impression that with time evolution, both the excita¬ 
tion energies and maximum fragments become small. In 
the case of central and semi-central collisions, the inul- 
tifragmention process can be obviously observed. The 
collective incident kinetic energy has transferred to the 
excitation energy due to the collisions between two nu¬ 
clei with same size. Then the excited fragments decay 
by emitting free nucleons and light charged particles. So 
the yield of fragment with charge Z < 10 increase with 
time evolution. In the case of peripheral collisions, af¬ 
ter t = 500 fm/c, there is only a few intermediate mass 
fragments (IMFs) yielded. The spectators with very high 
excitation energies are formed. With time evolution, the 
spectators decay by emitting free nucleons, and it only 
increases the yield of fragment with change around 60, 
but has a little contribution to the yield of IMF. Even 
until 3000 fm/c, there are still certain number of frag¬ 


ments with high excitation energies. It indicates that 
the statistical description is necessary for the peripheral 
collisions of this reaction. 

By using the hybrid model of ImQMD-l-GEMINI, we 
study the charge distribution of fragments in the reac¬ 
tions of ^®’^Au-|-^®^Au at incident energies of Ec.m. = 35 
MeV/u. The comparison between calculations and ex¬ 
perimental data are presented in Fig. 3. It should note 
that the incident energy is 40 MeV/u for experimental 
data of central collision. We do not think there should 
be quite difference to the data of 35 MeV/u. According 
to the experiment setup, in the case of central collision, 
products in the center of mass angles between 0° and 
90° are chosen to compare with the experimental data 
[31]. In the semi-central and peripheral collision, prod¬ 
ucts with charge up to the beam charge at 0iab from 3° 
to 23°, and products with charge up to Z=20 covering 
the angular range from 23° to 160° are chosen to com¬ 
pare with experimental data [30]. One can see from the 
figure that, in the case of central and semi-central col- 
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FIG. 2: (Color online) Time evolutions of the excitation energies per nucleon of primary fragments with charge number Z > 4 
for at incident energy of i?c.m. = 35 MeV/u. The average numbers of free nucleons at each event are also presented 

in each sub-figures. 


lision, the ImQMD model reproduces the experimental 
data reasonably well even without the statistical calcula¬ 
tions. In the early stage of reactions, there are still some 
large fragments with high excitation energies. With time 
evolution, the results directly from the ImQMD model 
approach gradually to those from the ImQMD-bGEMINI 
calculations, because the excitation energies of the frag¬ 
ments are exhausted by emitting particles. But in the 
case of peripheral collisions, the calculations without the 
statistical model obviously underestimate the yield of 
IMFs with 20 < Z < 40. It needs a much longer time 
for the heavy fragments to undergo large deformation and 
breakup. Because of the absence of some Fermionic prop¬ 
erties in the ImQMD model, some nucleons with high mo¬ 
mentum will be emitted more easily and take away some 
excitation energies. With smaller and smaller excitation 
energy, the heavy fragments become difficult to breakup. 
One can also see that, even the simulation stops at 3000 
fm/c, there is still obvious difference between the results 
with and without the statistical model being involved. It 


means that some fragments with high excitation energies 
will enter the secondary decay process. From the calcu¬ 
lations, we find that the experimental data can be well 
reproduced with the ImQMD-bGFMINI model at switch 
time t=500 fm/c. Finally what we want to mention is 
that, without any readjusting, our model can simultane¬ 
ously reproduce the data taken from two different exper¬ 
iments with very different apparatus; INDRA-I-ALADIN 
for the central collisions, and MULTICS for the periph¬ 
eral ones. It illustrates that our model has powerful pre¬ 
diction ability. We note that, in the present calculations, 
the angular momentum of fragments is not taken into 
account for the simplicity. In the fragmentation process, 
the emissions of nucleons and light particles could take 
away some angular momentum of heavy fragments. The 
neglect of angular momentum might be at the origin of 
the (small) discrepancy with experimental data which is 
still observed for the most peripheral collisions. The in¬ 
vestigate on the influence of angular momentum of frag¬ 
ments is under way. 
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FIG. 3: (Color online) Comparison between calculations and experimental data for charge distributions of fragments in 
^®^Au+^®^Au. The experimental data are taken from Ref. [31] for central collision, and Ref. [30] for semi-central and 
peripheral collision, respectively. 


IV. SUMMARY 

The multifragmentation in the intermediate heavy-ion 
collisions has been investigated by using the hybrid model 
of Improved Quantum Molecular Dynamics model plus 
the statistical model GEMINI. The excitation energy of 
fragments in heavy-ion fusion reactions is firstly tested 
with the ImQMD model, and the calculation results look 
quite reasonable. The roles of dynamical model and sta¬ 
tistical model in the description of the reaction mech¬ 
anism are investigated by analyzing on the fragments 
charge distribution of 35 MeV/u ^®^Au-|-^®^Au at various 
impact parameters. We find that, in the case of central 
and semi-central collisions, the calculations without the 
statistical calculations are in good agreement with the 
measured fragments charge distribution. The statistical 
process with nucleons and light charge particles emission 
can be partially represented by dynamical model. How¬ 
ever, in the case of peripheral collisions induced by heavy 


nuclei, the statistical model is still necessary to describe 
the multifragmentation. Otherwise the IMFs yields will 
be underestimated due to the absence of quantum cor¬ 
relations in the ImQMD model. By using the hybrid 
model of ImQMD-l-GEMINI, the charge distribution of 
fragments in Au-|-Au at 35 MeV/u can be reproduced 
reasonably well for both central and peripheral collisions. 
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